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Abstract 



We present results from CCD observations of comet 2P/Encke acquired at Steward Observatory's 
2.3m Bok Telescope on Kitt Peak. The observations were carried out in October 2002 when the 
comet was near aphelion. Rotational lightcurves in B, V and i?-filters were acquired over two nights 
of observations, and analysed to study the physical and color properties of the nucleus. The average 
apparent i?-filter magnitude across both nights corresponds to a mean effective radius of 3.95 ±0.06 
km, and this value is similar to that found for the V and B filters. Taking the observed brightness 
range, we obtain a/b > 1.44 ± 0.06 for the semi-axial ratio of Encke's nucleus. Applying the axial 
ratio to the i?-filter photometry gives nucleus semi-axes of [3.60 ± 0.09] x [5.20 ± 0.13] km, using the 
empirically-derived albedo and phase coefficient. No coma or tail was seen despite deep imaging of 
the comet, and flux limits from potential unresolved coma do not exceed a few percent of the total 
measured flux, for standard coma models. This is consistent with many other published data sets 
taken when the comet was near aphelion. Our data includes the first detailed time series multi- 
color measurements of a cometary nucleus, and significant color variations were seen on October 3, 
though not repeated on October 4. The average color indices across both nights are: (V — Ft) = 
0.39 ± 0.06 and (B — V) — 0.73 ± 0.06 (R = 19.76 ± 0.03). We analysed the .R-filter time-series 
photometry using the method of Harris et al. (1989) to constrain the rotation period of the comet's 
nucleus, and find that a period of ~11.45 hours will satisfy the data, however the errors bars are 
large. We have successfully linked our data with the September 2002 data from Fernandez et al. 
(2005) - taken just 2-3 weeks before the current data set - and we show that a rotation period of 
just over 11 hours does indeed work extrememly well for the combined data set. The resulting best 
fit period is 11.083 ± 0.003, consistent with the Fernandez et al. value. 
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1 Introduction 

Probing the physical properties of cometary nuclei is a challenge due to the presence of dust and gas 
comae that envelope the nuclei at small heliocentric distances. Cometary nuclei represent the least 
altered material from which to constrain the conditions and formation mechanisms of the early solar 
system as well as various subsequent evolutionary processes. By far the most well studied cometary 



population are the Jupiter-family comets (JFCs). This population are the most accessible due to 
their proximity to Earth and because many have small or negligible surface outgassing beyond 3 AU 
from the Sun, allowing the reflected and/or thermally emitted flux from the nucleus to be detected 
directly. The formation location of JFCs is most likely the trans-Neptunian region (Duncan et al. 
1995; Ip and Fernandez 1997). In contrast, the nuclei of long-period comets and Halley-type comets 
most likely formed closer to the Sun in the region of the giant planets. 

Studying the very low activity, near-dormant JFC comets provides a means for understanding 
the end states of comets. Indeed, comet Encke - a highly evolved member of the JFCs and the focus 
of this paper - may be close to becoming an inert comet. Thus it is important to study the surface 
and bulk physical properties of its nucleus. A systematic survey of such transitional objects will 
allow comparisons to be made with the so-called ACO population (asteroids in cometary orbits), 
or that of the near-Earth asteroids, to establish possible physical links. 

Full physical characterization, including determination of spin-rates has been accomplished for 
only a few comets. For a recent and thorough review of this topic see Lamy et al. (2004) as well 
as reviews by Samarasinha et al. (2004) and Weissman et al. (2004). Rotation statistics of comets 
are important for studying the physical evolution of these bodies and how they compare with their 
suspected progenitor population, the Kuiper belt objects. These distributions are expected to be 
different, especially due to outgassing torques that act on the nuclei as they enter the inner Solar 
system. Nucleus rotation and shape information can be obtained from time-series imaging of the 
comet, but unfortunately full nucleus lightcurves have been obtained for very few comets, due to 
the difficulties noted above. With full rotation lightcurves one also obtains a better measurement 
of the mean size of the nucleus, which aids in deriving a more accurate size distribution. 

Simultaneous multi-color lightcurves have not been obtained for any cometary nuclei and so the 
only data we have that points towards surface inhomogeniety comes from in-situ data acquired 
from spacecraft frybys, including Giotto, Vega, Deep Space 1, Stardust, and Deep Impact. When 
19P/Borrelly was imaged by the Deep Space 1 probe, the surface normal- reflectance was seen to 
vary from 0.01-0.05 (Buratti et al. 2004). 

Comet Encke has the shortest orbital period of any known comet (3.3 years), and its aphelion 
distance of 4.1 AU is well within the orbit of Jupiter. This orbital stability and short periodicity 
has allowed observations over many apparitions. Encke is therefore one of the most well studied 
comets, and is one of the few comets for which the albedo and phase coefficient have been measured 
(Fernandez et al. 2000). The nucleus has one of the steepest phase darkening slopes measured for 



a cometary nucleus with a phase coefficient of 0.06 magnitudes/ degree, measured over the phase 
angle range of — 106°. The visual geometric albedo is also relatively high at 0.05 ± 0.02, although 
the uncertainty overlaps the commonly accepted canonical albedo of 0.04. 

Earlier attempts were made at constraining Encke's rotation period. Time-series optical photom- 
etry from Jewitt and Meech (1987) show a most likely period of 22.43±0.08 hours, whereas Luu and 
Jewitt (1990) found a best-fit period of 15.08 ±0.08 hours. Both papers indicate that other periods 
were also consistent with the data. Thermal infrared time-series observations by Fernandez et al. 
(2000) are also consistent with the 15.08 hour period. A large data set was presented by Fernandez 
et al. (2005), acquired from July 2001 to September 2002, during which time the comet was near 
aphelion. They found that the synodic spin-rate is either 11.079 ± 0.009 hours or 22.158 ± 0.012, 
and that these periods were incompatible with the earlier estimates from Jewitt and Meech (1987) 
and Luu and Jewitt (1990), and vice versa. Belton et al. (2005) performed a detailed analysis of 
the available i?-band and 10/mi photometry and suggested that the nucleus may be in a complex 
or excited rotation state. Out of the two possible short-axis mode (SAM) states that provide good 
fits to the periodicities, Belton et al. found that the most likely state has a precessional frequency 
for the long axis about the total angular momentum vector P^ of 11.1 hours, and an oscillation 
frequency around the long axis P^ of 47.8 hours. 

The Encke nucleus was observed using the Arecibo radar during the close approach in November 
2003 and is thus the first comet to yield radar detections at multiple apparitions (Kamoun et al. 
1982; Harmon and Nolan 2005). The new radar data supports the recently reported 11 hour rotation 
period, and excludes the longer 15 and 22 hour periods that were previously suggested. Harmon 
and Nolan combined both radar and earlier IR data to obtain a solution for Encke's shape and size 
giving a mean effective radius of 2.42 km and an unusually large axial ratio of 2.6. Finally, comet 
Encke is also one of the very few comets to have its dust trail imaged at visual wavelengths (Lowry 
et al. 2003; Sarugaku et al. 2005). 

Here we present results from new CCD observations of comet Encke, acquired at Steward Ob- 
servatory's 2.3m Bok Telescope on Kitt Peak. The data includes the first detailed time series 
multi-color measurements of a cometary nucleus. The observations were carried out in October 
2002 when the comet was near aphelion at a heliocentric distance of 3.9 AU, and modest phase 
angle of ~ 7.2°. This epoch was close to some of the Fernandez et al. (2005) measurements, and 
so we attempt to link the two data sets to refine our spin-rate result. In section [2] we describe our 
observations in detail, as well as the lightcurve extraction techniques that were used. In section [3] 



we investigate the time-series color measurements. This is followed by a discussion of the derived 
size and projected shape of the Encke nucleus in section HI Section [5] includes a detailed discussion 
of the rotation properties that have been extracted from the photometry, and our efforts to link 
them with the Fernandez et al. measurements, taken close in time to our observations. A summary 
of our main conclusions is presented in section [61 

2 Observations and lightcurve extraction 

CCD imaging of comet 2P/Encke was obtained using Steward Observatory's 2.3m Bok telescope on 
Kitt Peak, on the nights of October 3 and 4, 2002 (UT). The images were obtained using an NSF 
Lick 3 2048x2048 pixel CCD at the telescope's Cassegrain focus. The pixel scale was 0.30"/pixel 
in 2x2 binned mode, and the total field of view was ~ 5.1 x 5.1 arcmin. We used the Harris B, 
V and -R-filters which have peak-transmission wavelengths at approximately 440 nm, 525 nm, and 
590 nm, respectively. These filters are roughly equivalent to the Johnson (B, V) and Kron-Cousins 
(R) filters. 

Our aim was to obtain time-series i?-filter photometry to determine the rotational properties of 
the nucleus from the lightcurve, and to obtain additional time-series imaging in the V and B filters 
in order to assess the colors of the object and to look for possible color variations as the object 
rotated. The observing strategy was to cycle through the three broadband filters by repeating 
the following sequence: R — V — R — B — R — V — R — B — R... while stopping occasionally for 
standard star observations when the object passed across or near background stars. This sequence 
gave approximately twice as many R-fAter images as V and B, and thus denser temporal/phase 
coverage of the lightcurve for the rotational analysis. 

The night-1 observations (October 3) include 18xi?-filter images with exposure times of 400 
seconds, SxK-filter images with exposure times of 400 seconds, and 7x.B-filter images with slightly 
longer exposure times of 500 seconds due to the expected drop in S/N at this wavelength. The night- 
2 observations (October 4) include 20xi?-filter images, 9xV-fi\teT images, and 8x5-filter images, 
all with exposure times of 400-500 seconds. For image processing purposes, a set of twilight sky 
exposures were taken through each filter in use, and used to flat-field the bias-subtracted images. 
All other instrumental artifacts such as cosmic rays and bad rows/columns were removed in the 
standard manner. 

For each wavelength pass-band, the cometary instrumental magnitudes were measured and 



compared with the brightnesses of several non-variable background stars to extract the rotation 
lightcurve. This was applied to the B, V and i?-filter data. As all cometary images were taken with 
the telescope tracking at the comet's apparent rate of motion, a slight trailing effect was introduced 
to the background stars. However, the affect on the S/N of the lightcurve data points was negligible. 
These relative instrumental magnitudes were calibrated using observations of the Landolt standard 
field PG2213-006 (Landolt, 1992), taken at a range of airmasses throughout each night. Observing 
conditions were photometric on both nights. 

Aperture photometry was performed on the comet and background stars through apertures large 
enough to take in at least 99% of the flux (i.e. 3xFWHM of background PSF), therefore seeing 
changes between images had no effect on the extracted lightcurve. The extracted B, V and R- 
filter lightcurves are shown graphically in Figured] for October 3 and 4. The BVR-filter apparent 
magnitudes and corresponding mid-exposure UT-Day are listed in Table [2] The mean B, V and 
i?-filter apparent magnitudes are listed in Table [3] for each of the two nights of observation. All 
image processing, photometry, and calibration were performed using the IRAF program (Tody 1986, 
1993). 

We searched for signs of a resolvable coma around the comet by measuring the surface brightness 
profile of the comet and comparing with the stellar-background brightness profile (or PSF). As all 
cometary images were taken with the telescope tracking at the comet's apparent rate of motion, 
the cometary profile was measured by aligning all i?-filter frames on a given night on the comet and 
co-adding them, and then measuring the azimuthally averaged brightness profile through a series 
of circular annulli. As the background stars were slightly trailed in each exposure, the brightness 
profile was extracted using the method of Lowry and Fitzsimmons (2005). This method provides 
the effective brightness profile which would have been obtained under sidereal tracking. The profiles 
were indistinguishable. In Figure [2] we show the i?-filter co-added images from each night, and one 
can see that there is no evidence of a tail, or even the dust trail that was previously detected at 
optical wavelengths (Lowry et al. 2003; Sarugaku et al. 2005). The comet appears as a sharp point 
source. 

3 Color variations 

The imaging sequence discussed above was performed so as to simplify, as much as possible, the 
interplotation between apparent magnitudes from the three filters for the color determination. With- 



out adequate interpolation between niters, the resulting color measurements will be confounded by 
the effects of varying nucleus cross-section as the object rotates. Each V-filter image was bracketed 
by two i?-filter images. Therefore, to obtain the (V-R) color index at the time of the V-filter image, 
we took the average of the two adjacent i?-filter apparent magnitudes and subtracted this from the 
V-filter apparent magnitude value. As for the (B-V), each I?-filter image was bracketed by two V 
and i?-filter images in the following way: V — R — B — R — V, therefore we took to the average of 
the two nearest V-filter apparent magnitudes and subtracted the average from the 5-filter apparent 
magnitude value. Table [2] lists the color values obtained. The average color indices for October 3 
are: (V-R) = 0.40 ± 0.06 and (B - V) = 0.70 ± 0.08. The average color indices for October 4 
are: (V — R) = 0.38 ± 0.05 and (B — V) = 0.76 ± 0.07. Encke's colors are slightly redder than the 
equivalent Solar values of (V — R) Q = 0.36, and (B — V)© = 0.67, but are bluer than most other 
comets as well as previous values for Encke (Luu and Jewitt 1990). This is illustrated in Figure [3] 
where we compare our colors with those comets for which BVR photometry exists. 

We searched for possible variations of color with rotation by plotting the BVR color indices 
versus UT-Day and the results are shown in Figure HI No significant color variations are seen for 
October 4, but on October 3 we see a clear and systematic dip in (V — R) color index, which is 
not seen for (B — V). This systematic variation is considered real at the 2a confidence level. Also, 
there was no sign of time- varying atmospheric extinction, from inspection of the brightnesses of the 
background stars. However, if we are looking at a bare nucleus and the 11.089 hour rotation period 
from section 15.11 is the correct one - and we strongly suspect that it is - then color variations should 
be repeated from night to night. 

We explore the possibility that a brief outburst in gaseous emission occurred around Oct. 3.2 UT. 
C2 molecules within cometary coma have prominent emission bands at wavelengths that overlap the 
passbands of the B and V broadband filters. Brief outbursts in activity were recently detected by 
the Deep Impact probe on approach to comet 9P/Tempel 1 (A'Hearn et al. 2005). These outburst 
were due to active areas coming into view of the sun as the comet rotated. For some of the larger 
outbursts, residual activity was seen for up to 18 hours after the time of peak brightness. In the 
case of Encke, such an outburst would need to be a singular, transient event, not repeated on 
subsequent rotations of the nucleus as the active spot again comes into sunlight, at least during 
the two nights of our observations. This would satisfy the criterion that if the ~11 hour period is 
correct then we should be looking at roughly the same hemisphere of the nucleus on October 3 and 
4 and thus the color variations would be the same from night to night for a bare nucleus. If these 



color variations were due to an impulsive release of gaseous species then as the color-change event 
lasted approximately 2.4 hrs, the gas would need to completely traverse our photometric aperture 
within this time. At a geocentric distance of 3.02 AU, the gas expansion velocities would need to 
be at least 1.2 km s _1 for our chosen aperture sizes. This is at least a factor of ~ 3—4 times larger 
than the expected expansion velocities for these species at a heliocentric distance of 3.93 AU, where 
the expected blackbody sub-Solar temperature on the Encke nucleus would be ~ 204 K. Therefore, 
this scenario is very unlikely. Also, the active-comet scenario is further weakened as no coma or tail 
was seen in our deep images of the comet, nor in previous near-aphelion imaging. 

We can conclude that the color variations are significant and most likely imply inhomogeneity on 
the surface of Encke's nucleus (or be due to a combination of nucleus-surface and coma properties). 
Whether there is coma present or not, we do demonstrate that color variations of this order are 
possible to detect. Complex rotation, in conjunction with a highly elongated nucleus with significant 
surface inhomogenieties, may well explain both the color variations and the erratic photometric 
behaviour of the comet when it is near aphelion (see section H]) . 

As a note for future comet-nucleus observers - we feel that it is important to focus on .R-filter 
photometry for determining the rotation period from time-series photometry, to avoid these emission 
bands as coma may be present. Had the V or B filter been the only filter used, then our period 
analysis may have been adversely affected. We encourage future lightcurve investigations to include 
full rotational phase coverage at multiple bandpasses, to search for possible color variations on 
cometary nucleus surfaces, and to give added dimension to the rotation period analysis. 

4 Nucleus size and shape 

We calculated the mean nucleus radius in each of three filters for both nights using, 



where m\ is the average B, V or R apparent magnitude for a given night, A\ is the geometric 
albedo in the respective filter, [m] is the nucleus effective radius, Rh [AU] and A [AU] are the 
heliocentric and geocentric distances, respectively, a and /3 are the phase angle and phase coefficient, 
respectively, and m & \ is the apparent B, V or R magnitude of the Sun. For these calculations we 
use m Q R = -27.26, (V — R) Q = 0.36, and (B — V) & = 0.67. The derived mean radii for each filter 



A x r 2 N = 2.238 x iq 22 /? 2 A 2 10 a4(m ©*- mA+a/3) 




and each night are listed in Table [31 We use the empirically derived values for the geometric albedo 
and phase coefficient of 0.047 and 0.06 magnitudes/ degree, respectively (Fernandez et al. 2000). 
The average i?-filter apparent magnitude across both nights is 19.76 ± 0.03, which corresponds to 
a mean effective radius of 3.95 ± 0.06 km. 

We can scale this average .R-filter apparent magnitude to the geometry of the September 2002 
photometry listed in Fernandez et al. (2005), using the standard 5log(RhA) scaling law for an inert 
body, along with the linear phase- darkening correction using the phase coefficient above. We find 
that the scaled October magnitude is 19.40±0.03 which is very close to the Fernandez et al. average 
i?-magnitude of 19.34 ±0.02 (or mean radius = 4.06 ±0.04 km), and agrees at the 2a level, implying 
an inactive nucleus on both datasets. A la agreement can be achieved with an increase in the phase 
coefficient of just 0.01. Alternatively, if one assumes that the nucleus is active and it has a typical 
dust-production dependence on heliocentric distance (i.e. the brightness would scale as 10log(Rh)), 
then the scaled October magnitude is at least 0.1 magnitudes too faint for agreement with the 
Fernandez et al. average -R-magnitude. Of course, the differences in heliocentric and geocentric 
distances between the two data sets are small and so scaling of this type cannot definitively rule out 
the presence of an unresolved coma. The mean-radius for the other data sets listed in Fernandez et 
al., taken much farther back in time than our October 2002 measurements, are as follows: July 19, 
2001 - 4.19 ± 0.05 km; August 10-13, 2001 - 3.24 ± 0.04 km; September 21-25, 2001 - 4.61 ± 0.05 
km; and October 6-8, 2001 - 2.36 ± 0.09 km. Again, these values assume a common albedo and 
phase coefficient. 

We now compare our mean i?-magnitude with other published data that include a good sampling 
of the i?-filter lightcurve. The photometry of Jewitt and Meech (1987) results in a mean-radius of 
3.75 ± 0.09 km, using their 1985 data set and using the empirically derived values for the albedo 
and phase coefficient. This value is close to ours of 3.95 ± 0.06 km. However, the time-series 
photometry from Luu and Jewitt (1990) imply a smaller mean radius of 3.28 ± 0.06 km. In three of 
the eight cases where .R-filter lightcurves were obtained - and for the radar±IR study from Nolan 
and Harmon (2005) noted in section CD- the mean radius was measured to be between 2.36 and 3.28 
km. For the remaining five cases, the mean radius falls in the range of 3.75 — 4.61 km. Although 
there is good agreement between most of the data sets, clearly the photometric behaviour of comet 
Encke is still perplexing, and may imply an erratically behaved unresolved coma with a very steep 
brightness profile (as the flux contribution from a shallow steady-state coma is only a few percent 
of the total measured flux), and/or complex rotation of a very elongated nucleus. There is no 



correlation between the observed mean brightness and the observing geometry or orbital position 
of the available data sets. This erratic behaviour, along with the color variations noted above, will 
require further data and complex modelling to fully understand (the latter being outside the scope 
of this particular paper). 

The observed magnitude range Am was 0.35 ± 0.05 and 0.40 ± 0.04 for October 3 and 4, respec- 
tively. This result uses the .R-filter data only as the lightcurve coverage was best for this wavelength 
passband. We can infer limits on the nucleus shape using: 

- > io°- 4Am (2) 
b 

where a, b are the semi-axes of the nucleus. This simple model assumes a bi-axial ellipsoid shape 
and uniform surface albedo, and ignores phase effects. Taking the maximum observed brightness 
range, i.e. from October 4, we get a/b > 1.44 ± 0.06 for the Encke nucleus. Applying this axial 
ratio to our .R-filter photometry gives nucleus semi-axes of [3.60 ± 0.09] x [5.20 ± 0.13] km. These 
dimensions, as with our measured mean radius above, could well be upper limits to the size, given 
the erratic behaviour of the photometry across the various published data sets, possibly implying 
outgassing. Although again, we point out that no resolved coma or tail has been seen in any of the 
above data sets, even with very deep-imaging techniques applied in most cases. 

5 Nucleus rotation rate 

We applied the method of Harris et al. (1989) to determine the rotation period from the time-series 
photometry. We focus on the .R-filter data to constrain the rotation period as we have roughly 
twice as many data points at this passband. This method involves fitting an nth-order Fourier 
series to the relative magnitudes, which is then repeated for a wide range of periods until the fit 
residuals are minimized. We created periodograms over the reasonable range of ~ 1-30 hours. This 
range includes all previously proposed periods for Encke noted in section [TJ When fitting model 
lightcurves to the data, we start off with low order fits to get a more accurate feel for where the 
prominent periodicities reside. We then increase the fit order to refine the dominant periodicity 
and its associated uncertainty, and also the quality of the fit. We stop increasing the order once the 
quality of the fit no longer improves. The dominant periodicity must be consistent with our initial 
visual inspections. 



This method was applied to the i?-filter data and the resultant periodogram from 2nd order 
fits is shown in Figure [5] (a). There are four main features at 7.74, 11.45, 16.30 and 22.90 hours, 
and the phased i?-band light-curves look believable at all four periods [see Figures [5] (b)-(e)]. The 
11.45 ±0.11 hour period produces the best fit, but this value is longer than previous determinations 
(Fernandez et al. 2005; Lowry et al. 2003). The measurement uncertainty is quite large, which is 
most likely due to reduced temporal coverage. The period at 7.74 hours does not satisfy the B and 
^-filter data, particularly V, given the very poor overlap of data points between the two nights. 
Also, there is the obvious fact that it is a single-peaked lightcurve, which is physically unrealistic. 
There is similarly poor overlap in the phased B and V data points for the 11.45 hour period (and 
22.90 = 2x11.45 hours), although to a lesser extent. These are attributable to the color differences 
between the two nights, noted in section [31 

Finally, the wide periodogram feature at 16.30 hours requires some discussion. We essentially 
rule out this period as it is most likely due to the period-fitting software slotting the two lightcurves 
end-to-end, leading to a false periodicity. Of course, this does not mean that such periodicities are 
not viable solutions, but in this case there are no prominent harmonic features associated with this 
period. Notice how wide the feature is. We speculate that this may be a merging of two separate 
features at 16.30 hours and perhaps one at 15.48 hours (i.e. 2 x 7.74 hours), where the curve seems 
to level off slightly. Even though this is close to the previously published period at 15.08 hours 
(Luu and Jewitt 1990), and a could explain our observed color differences between Oct 3 and 4, we 
conclude that both periods at 15.48 and 16.30 hours are likely pathological cases that appear to fit 
only because there is no overlap in the phased data from the two nights. Also, the persistence of 
the 11 hour period producing the best fit over most of the available data cannot be taken lightly as 
shown by Fernandez et al. (2005), and as we go on to show now when we link our data with some 
of the Fernandez et al. data set. 

5.1 Combining our data with Sep 2002 

The Fernandez et al. (2005) results are based on a very large data set of i?-filter magnitudes 
spanning July 2001 to September 2002. We consider here the data from September 10-16, 2002 
only as this was obtained close to our current data set taken just 2-3 weeks later. The Fernandez et 
al. September 2002 data were taken using the University of Hawaii's 2.2m telescope on Mauna Kea, 
and Steward Observatory's 2.3m Bok telescope at Kitt Peak, and contains 296 xR filter data points. 
The comet was at heliocentric and geocentric distances of 3.97 AU and 2.97 AU, and the phase angle 



was ~ 1.6°. Of course, if we combine the two data sets then the resultant period will be naturally 
biased towards the Fernandez et al. value as most of the data are from that set. The main aim of 
this execise is to refine our period from section [5] by combining our data with their September 2002 
data, i.e. taken approximately two to three weeks earlier than ours, thus substantially increasing 
the time-base, and verifying that a period just over 11 hours really can satisfy both data sets. The 
derived period must successfully phase the October data onto the September data (i.e. the phased 
brightness maxima and minima must coincide). 

The resultant periodogram from 2nd order fits is shown in Figure [6] (a). There are prominent, 
well defined features at 7.689, 11.083, and 15.153 hours, with a double feature at 22.219 and 22.708 
hours. As expected, there is no hint of a period near 16.30 hours (see section EJ). Plots (b-f) are the 
September and October data sets folded to each of the prominent periodicities. We have left the 
magnitudes as apparent magnitudes to separate the two data sets for clarity and to highlight the 
change in mean apparent magnitude between the two dates. This brightness change is consistent 
with an inactive nucleus on both occassions, but a tightly bound coma cannot be ruled out (see 
section Hj) . 

One can see from Figure [6] (b) that a rotation period of 7.689 hours phases the October data set 
reasonably well but does not work for the September data. In both cases the result is a physically 
unrealistic single-peaked lightcurve. The phased data for the 15.153 hour period is somewhat 
interesting in that although the folded data for September is quite poor, there is a hint of similarity 
in lightcurve shape between the two data sets (i.e. one peak followed by a slightly larger peak). A 
rotation period of 11.089 ± 0.003 hours is clearly the most likely period from the periodogram in 
Figure [6] (a). The scatter in the September data is reduced significantly and the October brightness 
maxima and minima line up extremely well. This period agrees very well with the Fernandez et al. 
value of 11.079 ± 0.009, although now the synodic period is known to a slightly higher degree of 
accuracy. 

This raises an interesting point about our period determination in section[5]of 11.45 ±0.11 hours, 
as well as determining rotation periods in general from data taken over a few nights. One can see 
that the 11.45 hour value is somewhat larger that the refined value of 11.089±0.003 hours. Although 
the two values are still consistent at the 3.3<r level, we feel that the difference is still sufficiently 
large to make the point that one must proceed with caution when extracting rotation information 
from lightcurves based on only 2-3 nights, and time-bases on the order of weeks may be needed 
to get a solid synodic period value. In this case the true synodic period was at the extreme of the 



uncertainty range of our 2- day data set. 

Fernandez et al. discussed the possibility of a period around 22.158 hours but that period doesn't 
show up in the periodogram of the combined data set. The nearest feature to this period is the 
double feature at 22.219 and 22.708 hours. This double feature appears as a result of the reduced 
rotational-phase coverage of the October data (from just two nights data as opposed to seven nights 
in the September data set). Almost-equally good fits can be achieved by simply shifting the phased 
October data left or right to line up with either of the two peaks in the September phase plot. This 
is more clearly seen in figures E^e) and (f). The September phasing with P=22.708 hours is not 
good, but is slightly better for P=22.219 hours. In neither case do the peaks line up as well as for 
the 11.089 hour value nor are the fit residuals as low. We therefore conclude that the best solution 
for the combined data set is 11.089 ± 0.003 hours. 

5.2 Comparison with published rotation periods 

As noted in previous sections, there are multiple lightcurves that have been published, including 
the large data set of Fernandez et al. (2005). In section [5j where we consider just our October 
2002 data, we have arguments that support a rotation period at 11.45 ± 0.11 hours. This period 
produces the best possible fit to the i?-filter data but results in an odd-looking lightcurve that is 
difficult to explain in terms of shape, and does not work well if the lightcurves from all three filters 
are considered collectively. This period disagrees with the Luu and Jewitt (1990) best-fit period 
of 15.08 ± 0.08 hours and the thermal infrared observations by Fernandez et al. (2000) which are 
consistent with the Luu and Jewitt period. Jewitt and Meech (1987) quote a most likely period of 
22.43 ± 0.08 hours. However, our result agrees well with the more recent Fernandez et al. (2005) 
data sets which place the spin rate at either 11.079 ± 0.009 or 22.158 ± 0.012 hours. In Figure [7] 
we phase our October data to each of the previous periods, including the 11.014 hour period from 
Lowry et al. (2003), which is a subset of the Fernandez et al. (2005) data. 

Despite the obvious differences in the best-fit periods across the various data sets, there are 
glaring commonalities in that all the periodograms show prominent features in the ranges ~7-8, 
~11-12, ~15-16 and ~22-23 hours. The actual rotation rate must fall in one of these ranges, and 
we conclude that the period just over 11 hours in the correct one, at least for the most recent data 
sets, and assuming simple relaxed principle-axis rotation. Not only is this the best-fit period for 
the September 2002 data set (Fernandez et al. 2005) and our October 2002 data set, but when 



they are linked together, it clearly offers the best solution. The period at 15.153 hours for the 
combined data set (Figure [6]) does not phase well at all. The only negative argument regarding the 
~11 hour period is the persistent misalignment of the brightness peaks when the data are phased 
to this period. An important point is that new radar data (Harmon and Nolan 2005) supports 
the ~11 hour rotation period, and excludes the longer 15 and 22 hour periods. Invoking complex 
rotation, as suggested by Belton et al. (2005) (see section CP, in conjunction with a highly elongated 
nucleus with significant surface inhomogenieties, may well be the only way to solve the spin-state 
problem of Encke. However, without futher data that would enable the necessary modelling to take 
place, we choose the single-axis-rotation approach for now and arrive at a best-fit synodic period 
of 11.083 ±0.003 hours. 



6 Summary and main conclusions 

We present results from simultaneous B, V and i?-filter photometry obtained on October 3-4, 2002 
at the Steward Observatory 2.3m Bok telescope on Kitt Peak. Rotational lightcurves in all three 
filters were extracted and analysed to study the physical and color properties of the nucleus. The 
main conclusions from this work are as follows. 

1. The average color indices were measured for both nights and the values are very similar. The 
values are (V - R) = 0.40 ± 0.06 and (B — V) = 0.70 ± 0.08 for October 3, 2002 (R = 
19.78 ± 0.03); and (V — R) — 0.38 ± 0.05 and (B - V) = 0.76 ± 0.03 for October 4, 2002 (R 
= 19.74 ± 0.03). Encke's colors are slightly redder than the equivalent Solar values, but are 
bluer than most other comets as well as previous values reported for Encke. 

2. On October 3 we see a clear and systematic dip in (V — R) color index, which is not seen for 
{B — V) on the same night, and is not repeated on October 4. This systematic variation is 
considered real at the 2a confidence level. However, for a nucleus rotation period of 11.089 
hours, color variations should be repeated from night to night. To explain the color variation, 
we explore the possibility that gaseous emissions from C2 molecules at B and V wavelengths 
occurred around Oct. 3.2 UT, but this is ruled out based on gas-expansion-velocity arguments. 
Also, no coma or tail was seen in our deep images of the comet, nor in previous near aphelion 
imaging. Thus, we tend to favor the no-coma scenario. Complex rotation, in conjunction with 
a highly elongated nucleus with significant surface inhomogenieties, may well explain both the 
color variations and the erratic photometric behaviour of the comet when it is near aphelion. 



3. The average apparent i?-band magnitude across both nights is 19.76 ±0.03, which corresponds 
to a mean effective radius of 3.95 ± 0.06 km. This value is similar to that found for the V- 
and 5-filter photometry. Taking the observed brightness range, we get a/b > 1.44 ± 0.06 
for the axial ratio of Encke's nucleus. Applying this axial ratio to the .R-filter photometry 
gives nucleus semi-axes of [3.60 ± 0.09] x [5.20 ± 0.13] km, for an inactive nucleus. These size 
measurements use the empirically derived values for the geometric albedo and phase coefficient 
of 0.047 and 0.06 magnitudes/degree, respectively (Fernandez et al. 2000). 

4. We analysed the .R-filter time-series photometry using the method of Harris et al. (1989) to 
constrain the rotation period of the comet's nucleus. We find that a period of 11.45±0.11 hours 
satisfies the data, however the errors bars are large. A rotation period somewhere between 
15-16 hours can work for the current data set in terms of the resulting lightcurve shape, which 
is reasonably consistent across the three observed wavelengths, and would explain the color- 
variation difference between each night. However, we believe that solutions in the 15-16 hour 
range are pathologic cases that brings the data from both nights together with no overlap. 

5. It is clear that a period just over 11 hours is the most likely value as this consistently produces 
by far the best fits (cf. Fernandez et al. 2005: Belton et al. 2005), and is supported by recent 
radar observations (Harmon and Nolan 2005). We have successfully linked our data with 
the September 2002 data set from Fernandez et al. (2005) - taken just 2-3 weeks before the 
current data set - and we show that a rotation period of just over 11 hours does indeed work 
very well for the combined data set. The resulting best-fit period is 11.083 ±0.003, consistent 
with the Fernandez et al. (2005) value. 

Acknowledgments 

We thank H. Campins and another anonymous referee for their very helpful reviews. This work 
was performed in part while the first author held a National Research Council Associateship Award, 
and in part at Queen's University Belfast. Also, this work was performed in part at the Jet Propul- 
sion Laboratory under contract with NASA. We acknowledge additional support from the NASA 
Planetary Astronomy Program. We thank Steward Observatory for granting time on the 2.3m Bok 
Telescope. IRAF is distributed by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. We acknowledge JPL's Horizons online ephemeris 
generator for providing the comet's position and rate of motion during the observations. 



References 

A'Hearn, M.F., Belton, M.J.S., Delamere, W.A., Kissel, J., Klaasen, K.P., McFadden, L.A., Meech, 
K.J., Melosh, H.J., Schultz, P.H., Sunshine, J.M., and 23 co-authors. 2005. Deep Impact: 
Excavating Comet Tempel I. Science 310, 258-264. 

Belton, M.J.S., Samarasinha, N.H., Fernandez, Y.R. and Meech, K.J. 2005. The excited spin state 
of Comet 2P/Encke. Icarus 175, 181-193. 

Buratti, B.J., Hicks, M.D., Soderblom, L.A., Britt, D., Oberst, J. and Hillier, J.K. 2004. Deep 
Space 1 photometry of the nucleus of Comet 19P/Borrelly. Icarus 167, 16-29. 

Duncan, M.J., Levison, H.F. and Budd, S.M. 1995. The Dynamical Structure of the Kuiper Belt. 
AJ 110, 3073-3081. 

Fernandez, Y.R., Lisse, CM., Ulrich Kaufl, H., Peschke, S.B., Weaver, H.A., A'Hearn, M.F., Lamy, 
P.P., Livengood, T.A. and Kostiuk, T. 2000. Physical Properties of the Nucleus of Comet 
2P/Encke. Icarus 147, 145-160. 

Fernandez, Y.R., Lowry, S.C., Weissman, P.R., Mueller, B.E.A., Samarasinha, N.H., Belton, M.J.S. 
and Meech, K.J. 2005. New near-aphelion light curves of comet 2P/Encke. Icarus 175, 
194-214. 

Harmon, J.K. and Nolan, M.C. 2005. Radar observations of Comet 2P/Encke during the 2003 
apparition. Icarus 176, 175-183. 

Harris, A.W., Young, J.W., Bowell, E., Martin, L.J., Millis, R.L., Poutanen, M., Scaltriti, F., Zap- 
pala, V., Schober, H.J., Debehogne, H., and Zeigler, K.W. 1989. Photoelectric observations 
of asteroids 3, 24, 60, 261, and 863. Icarus 77, 171-186. 

Ip, W.H. and Fernandez, J. A. 1997. On dynamical scattering of Kuiper Belt Objects in 2:3 resonance 
with Neptune into short-period comets. A&A 324, 778-784. 

Jewitt, D. and Meech, K. 1987. CCD photometry of Comet P/Encke. ,4J93, 1542-1548. 

Kamoun, P.G., Campbell, D.B., Ostro, S.J., Pettengill, G.H. and Shapiro, I.I. 1982. Comet Encke 
- Radar detection of nucleus. Science 216, 293-295. 

Lamy, P.L., Toth, L, Fernandez Y.R. and Weaver H.A. 2004. The sizes, shapes, albedos, and colors 
of cometary nuclei. In Comets II (M. Festou, H.U. Keller and H.A. Weaver, Editors), Univ. 



of Arizona Press, Tucson, pp. 223-264. 

Landolt, A.U. 1992. UBVRI photometric standard stars in the magnitude range 11.5-16.0 around 
the celestial equator. AJ 104, 340-371. 

Lowry, S.C. and Fitzsimmons, A. 2005. WHT observations of distant comets. MNRAS 358, 641- 
650. 

Lowry, S.C, Weissman, P.R., Sykes, M.V. and Reach, W.T. 2003. Observations of periodic comet 
2P/Encke: Physical properties of the nucleus and first visual- wavelength detection of its dust 
trail. LPSC 34, 2056. 

Luu, J. and Jewitt, D. 1990. The nucleus of Comet P/Encke Icarus 86, 69-81. 

Samarasinha, N.H., Mueller, B.E.A, Belton, M.J.S. and Jorda, L. 2004. Rotation of cometary 
nuclei. In Comets II (M. Festou, H.U. Keller and H.A. Weaver, Editors), Univ. of Arizona 
Press, Tucson, pp. 281-300. 

Sarugaku, Y., Ishiguro, M., Miura, N., Usui, F., and Ueno, M. 2005. Optical Observations of the 
Comet 2P/Encke Dust Trail. In Proceedings of Dust in Planetary Systems, LPI Contribution 
No. 1280, pl27. 

Tody, D. 1986. The IRAF data reduction and analysis system. In Proc. SPIE Instrumentation in 
Astronomy VI (D.L. Crawford, Eds.) 627, 733. 

Tody, D. 1993. IRAF in the Nineties. In Astronomical Data Analysis Software and Systems II, 
A.S.P. Conference Series (R.J. Hanisch, R.J.V. Brissenden, and J. Barnes, Eds.) 52, 173. 

Weissman, P.R., Asphaug, E. and Lowry, S.C. 2004. Structure and density of cometary nuclei. 
In Comets II (M. Festou, H.U. Keller and H.A. Weaver, Editors), Univ. of Arizona Press, 
Tucson, pp. 337-358. 



Table 1: Geometry of comet during observations 



UT Date 


R h [AU] 


A [AU] 


a [cleg.] 


Oct 3, 2002 
Oct 4, 2002 


3.925-3.925 
3.923-3.922 


3.020-3.021 
3.025-3.027 


7.09-7.16 
7.35-7.43 



Table 2: Observation log and B, V and i?-filter photometry 



UTt m x ± a Filter t (V - R) ± a (B-V)±a 

(midtime) [sec] 



October 3, 2002 
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± 


0.03 


R 
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0.53 


± 
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0.49 


± 
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± 
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± 


0.06 


B 


500 
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0.07 


0.68 


± 


0.08 
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± 


0.03 
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0.40 


± 
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± 


0.04 


V 
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0.41 


± 


0.06 


0.64 


± 


0.08 


3.1631 


19.77 


± 


0.04 


R 
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0.43 


± 


0.05 








3.1711 


20.77 


± 


0.05 


B 


500 


0.37 


± 


0.07 


0.62 


± 


0.07 


3.1872 


19.80 


± 


0.03 


R 


400 


0.30 


± 


0.05 








3.1946 


20.10 


± 


0.04 


V 
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0.29 


± 
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0.70 


± 


0.08 
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± 
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R 
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± 
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± 
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B 
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± 
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± 
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V 
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Table 2: Continued: Observation log and BVR photometry 



UTt m A ± a Filter t (V - R) ± a (B-V)±a 

(midtime) [sec] 
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f Light time corrected 

The relative magnitudes have been calibrated and placed on a standard magnitude scale (Landolt 1992). 



Table 3: Mean effective radii from BVR photometry 



UT Date Filter m\ ± a Radius ±<r [km] 



Oct 3, 2002 


R 


19.78 ± 0.03 


3.89 ± 0.06 


V 


20.19 ± 0.04 


3.80 ± 0.07 




B 


20.89 ± 0.05 


3.74 ± 0.09 


Oct 4, 2002 


R 


19.74 ± 0.03 


4.01 ± 0.05 


V 


20.14 ± 0.04 


3.93 ± 0.06 




B 


20.92 ± 0.05 


3.74 ± 0.08 
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Figure 1: Time series B, V and i?-filter photometry for October 3 and 4, 2002. 




Figure 2: Co- added i?-filter images of comet Encke from each night of observation, taken at Steward 
Observatory's 2.3m Bok Telescope on Kitt Peak. The comet is marked with a circle in panel (b). 
The comet was stellar in appearance on both nights and in all filters. 



Colors of cometary nuclei 



CO 




Figure 3: Comparison of our colors for comet 2P/Encke with those of other Jupiter-family nuclei 
for which BVR measurements were obtained together. Encke's colors are slightly redder than the 
equivalent Solar values, but are bluer than most other comets as well as previous values measured 
for Encke. Color data taken from review by Lamy et al. (2005). Solar colours are marked with a 
star. 
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Comet 2P/Encke - Oct. 4, 2002 
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Figure 4: Color indices versus time for each night of observation. No significant color variations are 
seen for night 2, but the beginning of night 1 shows a clear and systematic dip in (V — R) color, 
which is not seen for (B — V). This may be due to either imhomogenious nucleus surface properties 
or a brief outburst in gaseous emissions visible at B and V wavelengths around Oct. 3.2 UT. 
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Comet 2P/Encke - Oct. 3-4, 2002 [PERIODOGRAM] 

(a) 2nd order fits 
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Figure 5: (a) Periodogram from 2nd order fourier fits to the i?-filter data. The four most prominent 
features are marked, (b-f) The BVR magnitudes phased to each of these four periodicities, along 
with a phasing at 15.48 hours (see text). We find that the 11.45 hour period offers the best solution, 
but does not necessarily fit well across all three filters. When we combine our i?-filter data with the 
September 2002 data set from Fernandez et al. (2005), only a period near 11 hours offers a common 
solution (see section I5TTT) . Filled symbols - October 3rd; Open symbols - October 4th. i?-filter data 
are denoted by circles, V^-filter by triangles, and _B-filter by squares. 



2P/Encke - Sept 10-16 and Oct 3 and 4, 2002 



Comet 2P/Encke - Sep. 10-16 and Oct. 3-4, 2002 [PERIODOGRAM] 
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Figure 6: Our data from October 2002 was combined with i?-band data from Fernandez et al. (2005) 
spanning September 10-16, 2002, to improve the accuracy of the rotation period determination, (a) 
Several prominent features are seen in the periodogram, with the most prominent and well defined 
being the period at 11.083 ± 0.003 hours, (b-f) The October and September data folded to each of 
these five periodicities. The September data are shown as filled circles and the October data are 
shown as open circles. 
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Figure 7: Comparison of our photometry with published periods, (a) Lowry et al. (2003) - subset of 
Fernandez et al. (2005), (b) Fernandez et al. (2000), Luu and Jewitt (1990), (c) and (d) Fernandez 
et al. (2005). Symbol notation same as Figure [51 



